Anomalous High-Energy Spin Excitations in La2Cu04 
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Inelastic neutron scattering is used to investigate the collective magnetic excitations of the high- 
temperature superconductor parent antiferromagnet La2Cu04. We find that while the lower energy 
excitations are well described by spin-wave theory, including one- and two-magnon scattering pro- 
cesses, the high-energy spin waves are strongly damped near the (1/2,0) position in reciprocal space 
and merge into a momentum dependent continuum. This anomalous damping indicates the decay 
of spin waves into other excitations, possibly unbound spinon pairs. 



High-temperature superconductivity occurs when in- 
sulating layered antiferromagnets such as La2Cu04 and 
YBa2Cu306.i5 are doped with electrons or holes [ij. 
La2Cu04 and the other parent compounds of the high- 
temperature superconductors are Mott-insulating layered 
square lattice antiferromagnets with large exchange cou- 
plings. Doping these parent antiferromagnets results in 
a collapse of the antiferromagnetic (AF) order, a metal 
insulator transition and the formation of a superconduct- 
ing phase with a high transition temperature. Many be- 
lieve that the global phase diagram of the cuprates is 
explained in a strong coupling theory in which the spin 
correlations important for superconductivity are already 
present in the undoped parent compounds ^nS]- For ex- 
ample, Anderson [3] proposed that the insulating parent 
compounds are described by a "resonating- valence-bond" 
(RVB) model. In this model, the unpaired electrons on 
the Cu atoms pair with neighbours and are described by 
a superposition of singlet pairs. Hole doping introduces 
vacancies which can hop in a background of singlets. The 
pre-existing magnetic singlets may become the charged 
superconducting pairs in the superconductor [1, 2]. 

In this paper, we use inelastic neutron scattering (INS) 
to measure the collective magnetic excitations of the in- 
sulating parent compound La2Cu04 over a wide range of 
energy and momentum. Our experiments were performed 
using the MAPS spectrometer at the ISIS spallation neu- 
tron source of the Rutherford- Appleton Laboratory. The 
use of a spectrometer with a large angular detector cov- 
erage and crystals with good mosaic, has enabled con- 
siderable improvements on previous studies [4] . Thus we 
are able to observe new contributions to the magnetic re- 
sponse which carry significant spectral weight and high- 
light the strong fluctuations present in La2Cu04 at high 
energies. 

The collective magnetic excitations in La2Cu04 are 
due to correlations between the S = 1/2 unpaired spins 
associated with the Cu atoms. INS probes the spin-spin 
correlations as a function of wavevector q and energy hw. 
The neutron cross-section may be directly related to the 
imaginary part of the magnetic susceptibility x^\c[,uj)^ 
which is a measure of the strength of the excitations for 



a given (q, cc;). We specify wavevectors in terms of their 
positions in reciprocal space Q = ha.* + kh* + Ic'^ of 
the tetragonal (or square) lattice, where a ~ 3.8 and 
c ~ 13.1 A. Our sample consisted of five co-aligned sin- 
gle crystals of La2Cu04, with total mass of 47.5 g. The 
crystals were annealed in Argon at 1073 K for 48h. Sus- 
ceptibility measurements showed a Neel temperature of 
320 ± 5 K. INS data were collected at 10 K with incident 
energies of 160, 240, 450 and 600 meV and corrected us- 
ing the Cu^+ ^x^-y'^ form factor. 

Fig.[T]shows typical data in the form of constant energy 
slices plotted as a function of in-plane wavevector (/i, k). 
The high-intensity regions in Fig. [T] are produced when 
the spin- wave (SW) dispersion surface co'(q) (Fig. [IJa)) 
intersects with the measuring plane. Fig. [2] and Fig. [3] 
show energy- and wavevector-dependent cuts through the 
excitations which allow the spin-wave pole to be located. 
From these and other cuts, we are able to determine the 
SW dispersion throughout the Brillouin zone. This is 
shown in Fig. IH The results for the SW dispersion are 
similar to those obtained in previous INS and RIXS stud- 
ies |^,l5|]. 

An approach which has been widely used to describe 
the magnetic excitations in La2Cu04 is to assume that: 
(i) the magnetic interactions are described by a Heisen- 
berg model on a square lattice; (ii) the ground state has 
Neel order; and (iii) the excitations can be computed 
using classical (large- 6*) linear spin- wave theory (SWT). 
Following previous work [4] we used a Heisenberg Hamil- 
tonian with higher order coupling to model the magnetic 
interactions: 
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where J, J' and J" are the first-, second- and third- 
nearest-neighbour exchange constants and Jc is the ring 
exchange interaction coupling four spins. The exchange 
constants in Eq. [T] can be estimated in a Hubbard 
model in terms of the hopping energy t and double 
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FIG. 1: (Color online) (a) Dispersion as predicted by SWT 
(color indicates the scattering intensity), (b) The measured 
X''(q,cj) in units of fi% eV"^ f.u."^ at T = 10 K. At the 
highest energies (E>320 meV), the response is strongest near 
q = (1/2,0) and symmetry related positions. 



occupancy energy [3, H, Q- By fitting t and [/, 
we obtained J=143±2 meV, J'= J''=2.9±0.2 meV and 
Jc =58±4 meV. The resulting dispersion surface is shown 
in Fig.lIJa) and dispersion curves in Fig.HJ^a). A slightly 
improved fit to the intermediate-energy dispersion points 
can be obtained by including higher order hopping or ex- 
change terms as in Refs. [7, 8]. 

The SW dispersion has been measured in detail in 
a number of 5* = 1/2 square lattice Heisenberg anti- 
ferromagnets [1, The dispersion along the zone 

boundary between (3/4,1/4) and (1/2,0) is subject to 
two competing effects. For a nearest neighbor only cou- 
pling, SWT predicts there is no dispersion. Quantum 
fluctuations, not taken into account by SWT, cause the 
one-magnon energies to be renormalized raising (3/4,1/4) 
with respect to (1/2,0) This is observed experimen- 
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FIG. 2: (Color online) Magnetic excitation spectrum at vari- 
ous wavevectors in La2Cu04. (a) Ordering wavevector show- 
ing the two-magnon continuum, (b) (1/4,1/4) position on 
the magnetic zone boundary, (c) (1/2,0) position where the 
spin wave pole is strongly damped (zone boundary anomaly) . 
Lines show resolution-convolved fits to: a spin wave theory 
one- plus two-magnon cross-section [SWT(1+2M)]; quantum 
Monte Carlo simulation (QMC); and a generic power law con- 
tinuum lineshape (Eq. [2|). Data have been averaged over 
equivalent points for h,k < 1 and a background at (1,0) has 
been subtracted. 



tally when t/U (and J) is small 0-[l2|. For larger t/U 
(e.g. La2Cu04) the opposite happens as the higher 

order exchange coupling (Jc) in Eq. [T]raises (1/2,0) with 
respect to (3/4,1/4). 

SWT also makes a prediction about the q-dependence 
of the intensity of the SW pole. The solid line in Fig.Hl^c)- 
(d) shows this SWT prediction compared to the inte- 
grated intensity obtained from our data. Overall there 
is good agreement with the SWT prediction expect near 
q = (1/2,0) [or (7r,0)], where the SW peak intensity is 
strongly suppressed. Inspection of the energy-dependent 
scans Fig. [2fb)-(c) shows that in the region of q-space 
near (1/2,0), where the SW pole is suppressed, there is 
a significant high-energy tail to the magnetic response. 
This is also visible in the q-dependent plots Fig. [T] for 
£^>320 meV, which show the persistence of the scattering 
well localized near (1/2,0) and symmetry related points 
to the highest energies investigated. The high energy tail 
and intensity modulation can be seen most dramatically 
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FIG. 3: (Color online) Spin- wave and continuum (two- 
magnon) scattering in La2Cu04. (a)-(b) Cuts along Q = 
(/i, 0.5) for energy transfers /kj=232.5 and 280 meV and inte- 
gration ranges shown. The dotted line is a fit of SWT (one- 
magnon only) convolved the instrumental resolution with a 
quadratic background. The solid line is a fit of one-magnon 
plus two-magnon model. Right panels are magnified. 



in the image of the magnetic excitations along the high- 
symmetry directions shown in Fig. IH^b). Here, the tail 
at (1/2,0) and the intensity modulation between (1/2,0) 
and (3/4,1/4) can be seen. We also see evidence for con- 
tinuum scattering at lower energies inside the spin-wave 
cone as illustrated by the scans in Fig. [3l In an antiferro- 
magnet with collinear Neel order, a two-magnon contin- 
uum [nl, Q, [HI is expected due to the coherent creation 
of two independently-propagating spin waves with op- 
posite spin {Sz = ±1). This has been observed in other 
square-lattice antiferromagets such as Rb2MnF4 [161 and 
copper deuteroformate tetradeurate (CFTD) |11|. We 
have computed the two-magnon cross section based on 
the spin- wave dispersion determined by Eq. [T] using the 
standard expression 14, 15| (see solid lines in Figs. [2] and 
[3]). While the two-magnon simulation gives a good over- 
all description of the extra low-energy continuum inside 
the spin- wave cone, it fails to account for the high-energy 
lineshape. Inspection of Fig. [2fa)-(c) shows that the re- 
sponse at q = (1/4,1/4) is overestimated and that at 
q = (1/2,0) underestimated. 

We also compared our data with quantum Monte Carlo 
(QMC) simulations [13, [I^ 

which are shown as a dashed 
line in Fig. [2l The QMC simulation provides a good 
description of the q-dependence of the spin-wave pole 
amplitude, but does not fully account for the high en- 
ergy tail at q=( 1/2,0) in Fig. Efc). Finally, we fitted 
the complete lineshape (peak and tail) with the generic 



power-law decay function 



(2) 



where cjq is an onset energy (peak position), is 
a wavevector-dependent amplitude factor and is the 
Heaviside step function. This power-law function is a 
generalization of the continuum scattering lineshape [l8| 
of the ID Heisenberg AF chain. It gives a good pa- 
rameterization of the high-energy zone-boundary data, 
with T] = 0.93 ± 0.08 at q=(l/2,0), where the continuum 
scattering is strong, and a smaller r] = OAT ± 0.05 at 
q=(l/4,l/4), where the continuum tail is much weaker. 
Integrating Eq. [2] (for E < 450 meV) at the two posi- 
tions we obtain /(l/2,0)//(l/4, 1/4) = 0.90 ±0.11 show- 
ing that most of the loss in intensity from the SW pole 
is in the tail at (1/2,0). A reduced SW intensity and 
high energy tail at (1/2,0) have also been observed in the 
small t/U weak-exchange square-lattice antiferromagnets 
CFTD [nj and Cu(pz)2(C104)2 [13 • Although the tail 
appears to be weaker in CFTD than in La2Cu04 or 

Cu(pz)2(C104)2. 

In general terms, our results show that at the 
q=(l/2,0) position the spin waves are more strongly cou- 
pled to other excitations than at q=(l/4,l/4). This 
coupling provides a decay process and therefore damps 
the spin wave, reducing the peak height and producing 
the tail. The question is: "What are these other ex- 
citations?" An interesting possibility is that the con- 
tinuum is a manifestation of high-energy spinon quasi- 
Darticles proposed in theoretical models of the cuprates 



IHSI, Il3|, 11 91-12 Ij . These assume that Neel order co-exists 
with additional spin correlations with the magnetic state 
supporting both low-energy spin- wave fluctuations of the 
Neel order parameter as well as distinct high energy spin 
1/2 smnon excitations created above a finite energy gap 

2qI . 21 1 . Spinous are 5* = 1/2 quasiparticles which can 
move in a strongly fluctuating background. The anomaly 
we observe at (1/2,0) may be explained naturally in a 
model where spinous exist at high energies and have 
a d-wave dispersion [2o|, [HI with minima in energy at 
q=(±l/4, 1/4) and (1/4, ±1/4). Under these circum- 
stances, the lower boundary of the two-spinon continuum 
is lowest in energy at (1/2,0) and significantly higher at 
(1/4,1/4). This provides a mechanism for the spin waves 
at (1/2,0) to decay into spinous [with (1/4, ±1/4)] and 
those at (1/4,1/4) to be stable. 

The qualitatively new features in the collective mag- 
netic excitations observed in the present study are (i) 
a momentum dependent continuum and (ii) the q- 
dependence to the intensity of the spin- wave pole. We 
estimate the total observed moment squared (including 
the Bragg peak) is (M^)=1.9 ± 0.3 /i^. The continuum 
scattering accounts for about 29% of the observed inelas- 
tic response. The total moment sum rule [lij for 5* = 1/2 
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FIG. 4: (Color online) Wavevector dependence of the magnetic excitations in La2Cu04 showing anomalous response near 
(1/2,0). (a) One-magnon dispersion (T = 10 K) along lines in Fig. [JJc) (inset). Symbols indicate Ei: 160 meV (□), 240 
meV (A) and 450 meV (o). The solid line is a SWT fit based on Eq. [T] (b) Measured x'{^^^) throughout the Brillouin 
zone. Dashed circled highlights the anomalous scattering near (1/2,0). An energy-dependent background determined near 
(1,0) has been subtracted, (c) Experimentally-extracted one-magnon intensity. Line is a fit to SWT with renormalization 
factor Zcz=0.4±0.04. Dashed line in inset shows the AF Brillouin zone, (d) One-magnon intensity divided by SWT prediction 
highlighting anomalous loss in intensity at (1/2,0). 



implies (M^) = g'^ii\S{S + 1) = 3 We consider two 
reasons why we fail to observe the full fluctuating mo- 
ment of the Cu^+ ion. Firstly, our experiment is lim- 
ited in energy range to about 450 meV, thus there may 
be significant spectral weight outside the energy window 
of the present experiment. Raman scattering [22[ and 
optical absorption [i^ spectra show excitations up to 
about 750 meV. Recent resonant inelastic X-Ray scat- 
tering measurement also show high energy features [24] 
which appear to be magnetic in origin. The second rea- 
son why we may fail to see the full fluctuating moment 
may be covalency effects [25^, [1^. The Cu d^2_y2 and O 
Px orbitals hybridize to yield the Wannier orbital rele- 
vant to superexchange. This will lead to a reduction in 
the measured response. However, the (at most) 36% re- 
duction observed in La2Cu04 is substantially less than 
then the 60% reduction recently reported in the cuprate 
chain compound Sr2Cu03 |26|]. 

Our results have general implications for the cuprates. 
Firstly, they show that the collective magnetic excita- 



tions of the cuprate parent compounds cannot be fully 
described in terms of the simple SW excitations of a Neel 
ordered state. Secondly, they demonstrate the existence 
of considerable high-energy spectral weight above the SW 
pole. Thus, the ground state of La2Cu04 contains ad- 
ditional correlations not captured by the Neel/ SWT pic- 
ture which may take the form proposed by Anderson and 
others [sl, [l3|, ll9|, [20] . Deviations from SWT are strongest 
at high energies and are signalled by the observation of 
a q-dependent continuum of magnetic excitations and a 
strong q-dependent damping of the spin- waves. The SW 
pole is most heavily damped at the (1/2,0) [or (7r,0)] po- 
sition. This may be understood in a model where the 
spin waves are coupled to (and decay into) another type 
of dispersive excitation such as spinous. Angle resolved 
photoemission spectroscopy (ARPES) studies [13 show 
that doping La2Cu04 or other parent HTS compounds 
produces an underdoped superconducting state with a 
q-dependent d-wave pseudogap. The spinon dispersion 
that can explain the observed spin- wave damping in our 
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experiments has the same form as the pseudogap disper- 
sion, with a minimum at (1/4,1/4). 
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